Local properties in the two-dimensional t-t'-U model 
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We have studied the t-t'-U model by means ol the composite operator method. The effect of the 
bare diagonal hopping term t' that appears to be material dependent for high- Tc cuprate supercon- 
ductors is analyzed in detail. In particular, some local quantities are computed and a comprehensive 
comparison with the data by numerical simulations on finite size lattices is presented. The result 
show a good agreement with those obtained by Monte Carlo methods. 
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I. INTRODUCTION 

Since the discovery of the high-Tc superconductivity, 
there has been a great deal of discussion about the choice 
of an effective model suitable to describe the properties of 
the copper-oxide planes in the perovskite structure. Ex- 
tensive studies of the magnetic properties, showing one 
spin degree of freedom in the Cu-0 planeEl, have resulted 
in considerable evidence that the high-temperature su- 
perconductors may be modelled by an effective single- 
band model. In this line of thinking, one of the most 
studied model is the single-band Hubbard model which 
indeed can qualitatively describe many physical proper- 
ties experimentally observed in copper-oxide compounds. 
On the other hand, a particle-hole symmetric model can- 
not distinguish between electron- and hole-doped mate- 
rials. The addition of a finite diagonal hopping term 
t' has often been suggested to handle the complexity of 
the experimental situation for the cuprate being essential 
to reproduce various experimental observations. More- 
over, this electron- hole asymmetry in the next-nearest- 
neighbor hopping term, combined with a perfect symme- 
try of all the other effective parameters, emerges from 
various reduction procedures of multi-component elec- 
tronic models and seems to distingjiish the cuprates from 
a general charge-transfer insulatorH. 

In the next section the formulas for thad-t'-U model in 
the composite operator method (COM)0 framework are 
summarized. Our results and a comparison with dajta of 
numerical analysis by quantum Monte Carlo methoda arc 
also presented. Some conclusions are given in Sec. III. 



II. RESULTS 

Let us consider the t-t'-U model described by the 
Hamiltonian: 

H ^-li.'^c^ {i)c{i) -t^aijC^ {i)c{j) 



(1) 



where for a two-dimensional quadratic lattice with lattice 
constant a 



a(k) = i(cos(fc.a)-t-cos(M) 
(3 (k) — cos {kxtt) cos {kyu) . 



(2) 
(3) 



We use the spinor notation and drop the index of the 
spin freedom of electrons unless when it is necessary, 



(4) 



Following the COM ideas, we are interested in choos- 
ing a suitable asymptotic field for new bound states which 
appear, due to the strong correlations. Therefore, we in- 
troduce the doublet composite field operator 



with 



£.a (i) = Ca (i) (1 - («)) 



(5) 



(6) 
(7) 



The properties of the system are conveniently ex- 
pressed in terms of the two-point retarded thermal Green 
function: 



(8) 



In the static approximationEl, the Fourier transform of 
S{i,j) "is given by 



uj-m(k) J-i (k) 
where / (k) and m (k) are defined as 



/(k) 
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F.T. 



By considering a paramagnetic ground state, a 
straightforward calculation gives 
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mil (k) = -^/ii - 4i (A + a (k) (p 
-4t' (A' + /3(k)(p' + l-n)) 
mi2 (k) = TO21 (k) = 4t (A + a (k) (/; 
+4i'(A' + /3(k) (/22-p')) 
^22 (k) = (-/i + ?7) /22 - 4t (A + a (k)p) 
-4i'(A' + /3(k)p')- 
We use the following notation 

r (*) = ^«„v^(j-) 



A^(r (*))-(rr (*)r/t (^)) 

A' = (e^ (*))-(/ (*)^^ (0) 
1 
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(12) 

(13) 
(14) 
(15) 

(16) 
(17) 

(18) 
(19) 

(20) 



P = 4 ('^M («) '^M («)) - (c| (i) cx (i) (^cj {i) c\ (i) 
V=\ K {^) n, (*)) - L (^) c, {^) (c\ {^) c\ (z)) ''N (21) 



rip (i) = (i)(Tpc(i) being the charge (/x = 0) and spin 
(/.t = 1, 2, 3) density operator. 

The quantities A and A' are self-consistent parameters 
in the sense that they can be expressed in terms of the 
matrix elements related to the fermion propagator. The 
parameters p' and /i, the chemical potential, can be 
fixed by self-consistent equations 

n = 2(l-(e(z)^t(,))_<^(,)^t(,))) (22) 

(^(z)7?t(,)) = o. (23) 

The details will be presented elsewhere. The solution of 
the set of self-consistent equations allow us to compute 
the fermion Green function. 

We have computed the chemical potential and the dou- 
ble occupancy D — {n^ (i) (z)) for different values of 
the particle density, repulsive Coulomb interaction, tem- 
perature and bare diagonal hopping term t' . All the en- 
ergies are measured in units of t. In Figs.l and 2 our the- 
oretical results for n vs. ^ are presented and compared 
with the data obtainecLby a numerical study of a 8 x 8 
two-dimensional latticea. In Fig. 3 the double occupancy 
D is reported as function of the particle density. As it 
can be seen a negative t' decreases D when compared 
with t' = 0, while a positive value increases it. At half- 
filling the double occupancy is independent of the sign of 



t' as it is required by the symmetries of the model, and 
converges to the result for t' = 0. The agreement with 
the experimental data given in Ref. ^ is generally quite 
good. 



III. CONCLUSIONS 

By means of the COM, we have obtained a fully self- 
consistent solution for the t-t'-U model. As for the sim- 
ple Hubbard model, also in the case of the t-t'-U model 
our scheme of calculation can reproduce with good accu- 
racy the results of numerical simulation. In a forthcom- 
ing paper we shall continue the analysis of this model 
by considering the magnetic and transport properties 
which characterize the anomalous normal-state proper- 
ties of the cuprates. 
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